Methods for Estimating Capacities of Gaussian Quantum Channels 
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We present a perturbative approach to the problem of estimating capacities of Gaussian quantum 
channels. It relies on the expansion of the von Neumann entropy of Gaussian states as a function 
of the symplectic eigenvalues of the quadratures covariance matrices. We apply this method to the 
classical capacity of a lossy bosonic channel for both the memory and memoryless cases. 
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I. INTRODUCTION 

Quantum channels arc every means that convey quan- 
tum systems on whose states information is encoded. 
Formally it is a quantum map from input to output 
states [ij. The maximum rate at which information can 
be reliably transmitted through a quantum channel de- 
fines its capacity. Actually one can define several capac- 
ities depending on the kind of information transmitted 
(classical or quantum) and on the additional resources 
used in transmission [2j. 

Evaluation of quantum channels capacities is one of 
the most important and difficult problems of quantum 
information theory. Gaussian channels, which maps in- 
put Gaussian states into output Gaussian states, are 
among the simplest models allowing capacities investi- 
gation [U|. They are also relevant for experimental im- 
plementations [4[. 

A paradigmatic example of quantum Gaussian channel 
is the lossy bosonic channel where states lose energy en 
route from the sender to the receiver. The term bosonic 
arises because each input (respectively output) is repre- 
sented by a bosonic field mode. In turn, the effect of 
losses is usually modeled by letting each input mode to 
interact with an environment mode through a rotation 
(beam splitter) transform whose angle (transmissivity) 
determines the loss rate 

The classical capacity and the classical assisted capac- 
ity for such a channel were evaluated in Refs. 0, [|| by 
assuming each environment mode in the vacuum state. 
In such a case Gaussian states turn out to be the op- 
timal inputs. Then it has been proved that Gaussian 
inputs suffices to determine the quantum capacity 
More recently, it has been proven that they are optimal 
to reach the classical capacity when environment modes 
are in thermal states (including the vacuum) |8|. When 
more general states of environment are taken into ac- 
count, even non factorized ones which give rise to mem- 
ory effect [!|, the evaluation of capacities becomes much 
more demanding. 
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Here we develop a perturbative approach by expanding 
the von Neumann entropy of Gaussian quantum state as 
a function of the symplectic eigenvalues of its quadratures 
covariance matrix. The method turns out to be useful for 
characterizing capacities of Gaussian quantum channels. 
For the sake of presentation, we apply it to the classical 
capacity of a lossy bosonic channel. We then obtain lower 
bounds for the classical capacity, generalizing the results 
of Refs. 0,0. 

The paper is organized as follows. In Sec. II the main 
features of the lossy bosonic channel model are intro- 
duced. In Sec. Ill the general algorithm to find lower 
bounds on the classical capacity is presented. In Sec. IV 
the lower bound for the memoryless lossy bosonic chan- 
nel is given. In Sec. V the lower bound for the classical 
capacity is found for a memory lossy bosonic channel. 
Sec. VI is for conclusions. 



II. THE LOSSY BOSONIC CHANNEL 

Let us start from the consideration that in the evalu- 
ation of capacities of a Gaussian quantum channel one 
deals with the von Neumann entropy of Gaussian states, 
involving the function 
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where v is a symplectic eigenvalue of the covariance ma- 
trix characterizing a Gaussian state. Despite the absence 
of a small parameter for the expansion of the function g, 
we proceed as follows. The function g(v — 1/2) is not 
analytic in neighborhood of 1/2 and infinity, however af- 
ter subtracting the logarithm part it becomes analytic on 
the region v ^ 1/2. Hence, we may write (see also flol ]) 
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Then, to the zeroth-order approximation we have 

log 2 ,+ mT 
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where we have neglected terms of the order 0(l/v 2 ). Al- 
lowing perturbation of logarithm by the first terms in 
the series ([2]) we can also construct next-order approxi- 
mations. 

Below we restrict our attention to the classical capacity 
of a lossy bosonic channel. Since we shall consider Gaus- 
sian states, the action of this channel can be regarded 
as a linear map on the covariance matrices character- 
izing such states. We label covariance matrices acting 
on canonical quadratures (gi, . . . , q n ,pi, . ■ ■ ,p n ) by V"i n , 
V cnv and V out for input, environment and output states 
of the channel respectively. V c i/2 denotes the covariance 
matrix for the distribution of coherent amplitudes used 
to encode classical symbols. V" ut is the state of channel 
output averaged over encoding distribution (with covari- 
ance V c i/2). Also, we will refer to the eigenvalues of ma- 
trices Vi n , V cnv , Kut, V c \ and V ou t as 
and a u k respectively. Since n G N represents the num- 
ber of bosonic modes (channel uses), the dimension of all 
above matrices is 2n x 2n and k = 1, . . . ,n, while u de- 
notes either quadrature q or p (if u = q, then = p, and 
vice versa). In particular the lossy bosonic channel acting 
as a rotation (beam splitter) transform on the canonical 
quadratures, gives rise to the following relations among 
covariance matrices fill : 



Vout = V Mn + (1 - V) Venv, 
F out = V (Vin + V c l) + (1 - V) V* 



(4) 
(5) 



where rj represents the transmissivity of the channel. As 
usual (see e.g. [9J) we shall constrain the input energy by 
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2 -1r(V in + V tL ) = N+- i (6) 

with TV the average number of input photons per mode 
(channel use). 

An upper bound on the classical capacity is provided 
by the maximum of the Holevo information regularized 
in the limit of infinite channel uses 

C := lim C n ; C n := -maxx„ (7) 

n — >oo ft 

where the maximum is taken over all input states satis- 
fying the energy constraint © and is the Holevo-x 
quantity [H, ■ By coding theorems [H, , the quan- 
tity C turns out to exactly be the capacity for memoryless 
and forgetful channels [lj, [0 . 

For any n, one can look at n uses of the channel de- 
scribed by Eqs. ((U) and @ as a single channel acting 
on n modes. Then, a lower bound on the classical capac- 
ity can be computed by maximizing its one-shot capacity 
over Gaussian input states |16| 
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The Holevo-x quantity for the set G of Gaussian states 
reads [|[ 
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where Vf~ and v k are the symplectic eigenvalues of V out 
and Vout respectively. Taking the limit over n, we can as 
well define the quantity 



C := lim C r 

n — >oo 



(10) 



We first notice that by using fl3|) we have, at lowest 
order, 



k—1 



(ii) 



Then, we restrict our environment model to matrices of 
the form 
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with commuting blocks V qq and V pp . In this case we con- 
jecture that the maximum of x _ q uan tity ([9]) is achieved 
with matrices Vi n and V c i of the same form as (fl2|) . i.e. 
with null off-diagonal blocks. Furthermore, all diagonal 
blocks of all matrices are mutually commuting. This con- 
jecture is supported by numerical investigations relying 
on environment models of the form p2[) . As a conse- 
quence, we get the following expressions for symplectic 
eigenvalues 

(13) 



Vk — ^OgkOpk, Vk — ^/O-qkO-pk- 

Notice that both energy constraint ([6]) and symplectic 
spectrum f)13|) are preserved under orthogonal transfor- 
mations. Thus, without affecting the final result, we can 
consider all the involved matrices to be diagonal (see also 
the discussion in the appendix of [lq). 

Hence, the problem we are going to consider is the 
maximization ([5]) of the Holevo function © over Gaus- 
sian inputs under the constraints 

1 " 1 

y ^ [ink + iu*k + Cuk + Cu*k] = N + -, (14) 
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iuk > 0, 
Cuk > 0, 
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with Vk, Vk given by Eqs. (fT3l) . where 

o u k = Wuk + (1 - v) e uk, 

a u k = T\{iuk + Cuk) + (1 - r\)e uk . 



Here, the condition (fMj) is the input energy con- 
straint ©, while inequalities (|15p and (jTTJ) express the 
Heisenberg uncertainty relation. 



III. THE KARUSH-KUHN-TUCKER 
EQUATIONS 

The above maximization problem with constraints can 
be solved by considering eigenvalues of our matrices as 
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variables for Karush-Kuhn- Tucker (KKT) conditions in 
analogy with the same problem for classical Gaussian 
channels with memory [ITT ] . To simplify the calculations 
we will apply a technique similar to the simplex method 
in optimization theory. Assuming monotonic dependence 
on N of the eigenvalues {c u k} maximizing \ (which is 
supported by numerical investigation), we do not specify 
their sign in KKT-equations. At first, we will consider 
the formal solution of KKT-equations assuming negative 
eigenvalues as result, then we will specify how to elimi- 
nate them. 



Definition 1. Let us consider the eigenvalues c u k, 
u S {l,p} maximizing x as monotonically dependent 
functions of N £ [0,oo) ; then we will classify the k-th 
mode as belonging to the: 

• first stage iff c qk ,c pk = 0; 

• second stage iff c uk = 0, c Utrk > 0; 

• third stage iff c qk ,c pk > 0. 

By increasing N starting from zero we allow each mode 
to pass all stages in sequence from the first one up to the 
third one. 

The Lagrange function for n-modes channel reads 
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+ A y~ 22 [igk + ipk + £qk + Cpfe] - 2N -lj . 

Solving the KKT-conditions we find that the input state 
Vi n maximizing \ must be pure. When all modes belong 
to the third stage, the eigenvalues of Vo U t are all equal 
and the solution reads 



1 1 / e u k l — ?/ (TrV cnv 
c u k = N + - - - . / 1 — e u k I , 
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<7[(l-»7)(Ve^-l/2)], (20) 



where 
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is the average number of photons in the environment. 
The analytical lower bound in (f2"0"l) generalizes the ex- 
pression presented in [l6| . 



If we consider the lower bound ([20]) as a function on the 
set of memory models characterized by different V env each 
of them with the same M onv , we argue from Eq. (j2"0)) that 
the maximum is achieved when the state V env is pure. 

We now move to the general case. At first, we assume 
that the correct distribution of modes over stages is al- 
ready found. Then, suppose to know that the eigenvalue 
(f° r the mode h) belongs to the second stage and 
c u h = 0, it follows the transcendent equation for i u h 
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Vh ?; j (22 j 



dv h 



where c u *h = {a Ut h - (1 - v) e u,h)/i] - l/(4«„/,) and 
iu k h — l/(4? u 7 l ). This equation has real roots if both 
V h ,v h > 1/2, giving 

If 1 
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In turn, the requirement > Vh leads to the inequality 

V 



luh > 



^{a u „h - (1 - v) e u„h)' 



(24) 



Depending on the value of a Ut hi equation (|22[) can admit 
one root satisfying (|2"4"1) or none. 

Equation ([2^)1 can be formally written as the depen- 
dence a Ui h = fh{i u h)- Let A be the only parameter link- 
ing the KKT-equations of different modes, we define a 
new variable x 
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(25) 



getting a chain of equalities linking all modes of the sec- 
ond and third stages. Here modes m and I belong to 
the third stage, while modes h and t to the second stage 
(c p h = c qt = 0). Modes of the first stage are not included 
in P5|) and all give Vin-eigenvalues equal to 1/2. If some 
mode belongs to the third stage, its V; n -eigenvalues can 
be found from the relation (TIT))) . 

Taking into account stages discrimination, equation (J5]) 
can be rewritten as 



77Tr {2 , 3 | Cufc#0} y i; 
= [2n3 + U2]x, 



(1 - f?)Tr{ 2) 3|c u »540}Ku 



(26) 



where V- v 



Vin + V c \. Furthermore, n,- is the number 



of modes belonging to j-th stage (j = 1,2,3; n = m + 
77 2 +77,3) and Tr/2,3| c „n.^o} stands for the summation over 
all eigenvalues of second and third stages, except for the 
uk-th ones corresponding to c uk — 0. Also, the energy 
constraint © can be rewritten as 



Tr {2,3| Cufc #0}^in = 277 



N ■ 
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ink, (27) 
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where i uk — /jT 1 ^) and the double prime sum extends 
over ufc-th eigenvalues of the second stage, such that 
c uk = 0. Substituting Eq. (|2Tf into Eq. (|26|) we get a 
transcendent equation for the single variable x. Since 
all unknown eigenvalues can be expressed through x (see 
Eqs. (|25|) ) we can formally arrive at C_. 

To the zeroth-order approximation ([3]) the relation 

iuk = fk X { x ) ( see E( l- <E2>) S ives 

.(o) \/r] 2 + 16(l-r])xe Utk -r] 
t ° k ~ t <* = 8(1 ' (28) 

which allows us to express C_ as function of solution of 
only one algebraic equation for one variable x. By consid- 
ering the term l/v 2 of the decomposition ([2]) in Eq. ([22]) . 
we obtain the first-order approximation for the relation 
iuk = fk~ ( x )- Since Eq. (|22p cannot be exactly solved 
within this approximation, we solve it in neighborhood 
of the zeroth-order solution as linear perturbation. Thus, 
substituting + e uk in Eq. (|22|) instead of i uk and 
solving for e uk we obtain the first-order approximation 

ink =iul+ e uk, where 

£uk = {x — o u k)(x - o u ^k)ouki^l 

x{[2( 2 fe + x 2 - (o uk + o Uik )x) + (12(4 + 1)4]^ 
-2(l + 12^ 2 K fe .x}- 1 (29) 

where all eigenvalues are calculated to zeroth-order ap- 
proximation. 

Finally, let us discuss on how to find the correct dis- 
tribution of modes over the stages. As far as KKT- 
conditions themselves do not provide effective method 
to find stages distribution, we need to use some a 
proiri properties to write an effective algorithm. Since 
we have conjectured (supported by numerics) that c uk - 
eigenvalues are monotonic functions of N, the optimum 
for each mode has to be the third stage, or the second 
if the third is not admissible, or even the first if also the 
second one is not admissible. The third stage implies the 
existence of a solution inside some volume, the second 
stage a solution on the surface of that volume, while the 
first stage a solution on the edges of that surface. If the 
value of N is large enough, the solution is always inside 
a volume, so that all modes are in the third stage. By 
decreasing the value of N, the eigenvalues c uk pass from 
the volume to the wall, and then to the edges. On the 
basis of these arguments we can develop two algorithms: 
static and dynamic. 

The first step for both algorithms is the same: we check 
the positivity of all c uk we have found through the rela- 
tions (US]). If all of them are positive, C_ is given by the 
explicit analytical relation ([20]) and the problem is al- 
ready solved. If this is not the case, we move on to the 
next steps. 

• The static algorithm gives the following continua- 
tion on the bases of c M fc-signs found in the first step. 
We have to ascribe the third stage to modes with all 



positive c uk -eigenvalues, the second stage to modes 
with only one negative eigenvalue c uk and the first 
stage to modes with both negative c uk . All neg- 
ative c uk should be marked as already found and 
equal to zero. According to this distribution the set 
of KKT-equations should be formally solved (once 
again). If some second stage modes do not have 
a solution according to Eq. (|2"2"]) they have to be 
marked as belonging to the first stage. Here we ne- 
glect the condition ([24| for second stage modes, as 
violation of it implies negative c uk , thus their stage 
will be changed to the first one by the algorithm 
itself. If all found c uk are finally non-negative the 
problem is solved. If it is not, the procedure is iter- 
ated up until all c uk will become non-negative. This 
algorithm can be then considered as a consecutive 
correction of stages distribution. 

• The dynamic algorithm continues after the first 
step as follows. We solve the transcendent equa- 
tion on x by choosing different stages distributions 
for different values of x (during iterations) . In par- 
ticular, each mode at beginning is calculated as be- 
longing to the third stage for the current value of x. 
If this leads to negative c uk , the mode is marked as 
belonging to the second or the first stage in anal- 
ogy with the static algorithm (if some second stage 
modes violate Eq. (|24p their stage should be marked 
as the first one). Thus, we make stages distribution 
admissible, in quantum sense, for every value of x. 
The distribution of stages corresponding to a root 
of the equation for x is considered to be valid and 
can be used to calculate the correct eigenvalues. 

Both static and dynamic algorithms always yield the 
same stages distribution. The difference between the two 
algorithms can be clarified as follows. We have two ef- 
fective unknown "variables" for KKT-equations: stages 
distribution and x. One of these variables has to be set 
as internal and the another one as external during max- 
imization of Xn- The static algorithm uses x as internal 
variable, while the dynamic algorithm uses stages dis- 
tribution for that. Since dynamic algorithm is usually 
faster, below we will make use of it. Also notice that the 
eigenvalues i uk are always obtained as positive through 
these algorithms, therefore we did not specified their pos- 
itivity in KKT-conditions. 

If number of channel uses tends to infinity the dis- 
cussed procedure can be properly generalized by chang- 
ing the transcendent equations to equations on functions 
(spectral densities). However, if the considered model has 
some simple symmetries over stages, the general solution 
can be simplified by considering some parameters which 
mark the boundaries of stages. Below, we will show an 
example along this line. 

The solution of KKT-conditions can be interpreted as 
"quantum water filling" in analogy with usual (classical) 
"water filling" introduced for classical Gaussian channels 
with memory (see e.g. [l7jV The dependence of the found 
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spectral densities (also symplectic ones) from N is simi- 
lar to filling a vessel with water. The form of the vessel is 
defined by the model V cnv and transmissivity r\. The sym- 
plectic spectral density 77^ goes always up by increasing N 
(with respect to v^(N = 0)), while goes always down 
(or does not change). For environment models showing 
correlation (memory) among modes, the presence of the 
second stage gives rise to capillary effects on the edges 
of the vessel resulting to a "water level" with meniscus 
form. 

Below we apply the developed approach to channels 
whose environment is described by a covariance matrix 



V mv = N, 



1 









-sO 



(30) 



where N cnv is the average number of thermal photons per 
mode in the environment, f2 is a real symmetric nxn ma- 
trix and s € R is a parameter describing the environment 
properties. 



IV. THE MEMORYLESS CHANNEL 

We first consider the case of fi = I describing a memo- 
ry less channel with a thermal and squeezed environment. 
The first stage for one channel use is always degenerate 
and corresponds to N — and C = 0. The case of the 
third stage, which includes values of N ^> s, results to 
an explicit formula for the lower bound on the classical 
capacity (see Eqs. ([T5 )) - (jl20")0 : 



C = g[r)N + (1 - r])M env ] - g[(l - n)N cnv ], 



where, from the definition (121 



\[,. ny = j N CDV + ^ j COshs - i. 



(31) 



(32) 



This analytical result generalizes that of 0, Q . 

Considering the lower bound (|3"Tj) as a function on the 
set of environment models with fixed M env , we see that 
it shows a symmetry breaking over quadratures. In fact 
despite the symmetry of all equations over quadratures, 
the maximum of lower bound is achieved when e q ^ e p . 

The solution for the second stage is c u = 0, c Uji = 
2N + 1 — i u — l/(4i u ), i Ut = l/(4i„), where i u is a root 
of the transcendent equation 



J_ _ _±_\ v dg{V-l/2) 
o u a Ui / dv 



1 



dg^-l/2) 
dv 



(33) 



The condition c Ut > restrict the admissible region for 
i u to the interval 

N + i - y/N 2 + N <i u <N+- + ^N 2 + N. (34) 



To zeroth-order approximation ([3]) the solution of 
Eq. (33]) is 



7 -(o) _ 



i r 



v/l + (2iV+ 1)0 + 2 /4 -0/2 , (35) 



where 4> = e u 1? 7/(l — r /)- Hence, the classical capacity 
with this approximation is expressible in an explicit form. 

By considering the term 1 /v 2 of decomposition Q in 
Eq. (|33p . we obtain the equation for the first-order ap- 
proximation. Proceeding in analogy with the general 
case (see Sec. Ill) we find the first-order approximation 
,,,, where 



_ ,-(o) 



;(0) / _ \ 



2[^( 2 



■<£ ol{l2v* + l)] 



(36) 

and all eigenvalues in Eq. (f3"6"]l are calculated through 
zeroth-order approximation. Notice that the first-order 
approximation considered in Sec. Ill does not coincide 
with iu^ found. Actually, we have to jointly solve two 
equations in the case of general method applied to one 
channel use: the transcendent equation (|22[) and the 
equation for x (see Eqs. (|26l) and ([27]) ). where only 
Eq. (|22[) is approximated. However, the equation for x 
becomes analytically solvable in the case of one channel 
use, allowing us to reduce the original two equations to 
a single one given by (|33j). 

The formal algorithm to find the lower bound on the 
capacity is similar to general case considered in Sec. III. 
At first we solve KKT-conditions by considering the 
mode to be in the third stage. If one of the found 
c M -eigenvalues is negative we mark it as already found 
and set equal to zero. Then, we solve again the KKT- 
conditions by considering the mode to be in the second 
stage. This always brings us to the correct eigenvalues 
allowing us to find the lower bound on the capacity. 

Writing the input covariance matrix in the same form 
of ([30)1 (with Q = I), by the replacements N cnv — > 7V in 
and s — > r, we can relate the optimal degree of input 
squeezing r opt to the degree of environment squeezing s. 
It follows from Eqs. (TT5]) and ([33)1 that r opt = s for the 
third stage and 



ropt w = sign(s) In |yi + (27V + 1)0 + <f> 2 /4 - 0/ 

for the second stage in the zeroth-order approximation. 
At the transition point between stages there is a kink in 
the function r op t(s) (see Fig[T]-right). It reflects the fact 
that different stages correspond to solution of different 
systems of equations. 

The lower bound C_ found by the exact and approxi- 
mate analytical solutions is shown in Fig[l|-left as func- 
tion of s and for different values of 77. One can see that 
in the limit of large s the lower bound on the capacity 
only depends on the energy constraint N, specifically (see 
also pit) 



lim C(s, r), N, N cm ) = log 2 (27V + 1), 



(37) 
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FIG. 1: Lower bound on the classical capacity C_ (left) and 
optimal input squeezing r op t (right) vs s, for values of r\ go- 
ing from 0.1 (bottom curve) to 0.9 (top curve) with step 0.2. 
The values of the other parameters are N = N cnv = 1. On 
the left, exact and approximate solutions almost coincide. On 
the right, for each value of rj, the zeroth-order approximation, 
the first-order approximation and the exact solution are plot- 
ted (bottom to top). Solid and dotted parts of the curves 
correspond to the third and second stages, respectively. 




1.5 1 1 1 1 ■ 

0.4 0.6 0.8 Tl 0.4 0.6 0.8 



FIG. 2: Lower bound on the classical capacity (Z( s iv onv ) (left) 
and optimal environment squeezing Sat o[iv (right) vs rj for val- 
ues of iVenv going from (left curve) to 2 (right curve) with 
step 0.5. The value of the other parameter is N = 1. 



explaining why all curves C_ v (s) flow together to the same 
value when s — + oo (see FigJTJleft). The optimal degree 
of squeezing in the environment s* (which maximizes the 
lower bound on the capacity) is shown in Fig|2]-right and 
FigHJ-right as function of r] for different values of N env 
and N respectively. The lower bound on the capacity 
corresponding to these values of s* is plotted in FigJ^l-left 
and Figf4]-left. The critical transmissivity rf , represent- 
ing a threshold between transmissivities with finite and 
infinite values of s*, is drawn in Fig[3]-right as function 
of iVenv for different values of N. 

The points (a q , a p ) and [p q , o p ) are plotted in FigJ5]-left 
versus N for fixed values of s, rj and N cnv . The result- 
ing curves shows the geometry of the stage transitions 
and visualize the "quantum water filling" effect for one 
channel use. 



V. A CHANNEL WITH MEMORY 

Let us now consider an environment model (|30|) with 

= Si,j+t + Sij-i] i,j = !,•••, n (38) 




FIG. 3: On the left, quantities (o q , o p ) and (a q , a p ) vs N. The 
values of other parameters are iV cnv = s = 1, r) = 0.6. The 
second stage is marked by dotted curves. The first stage is a 
single point at N — 0. The whole third stage for the quantity 
(o q , o p ) is mapped into a single point, as Vi n does not depend 
on N in the third stage. On the right, critical transmissivity 
rf vs iVcnv for values of iV = 1, 500, 1000 from bottom to top. 
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FIG. 4: Lower bound on the classical capacity C_(s* N ) (left) 
and squeezing parameter s% (right) vs rj for values of N going 
from 0.5 (bottom curve) to 2.5 (top curve) with step 0.5. The 
value of the other parameter is iV cnv = 1. 



describing a lossy bosonic channel with memory 
Here s represents the degree of correlation among en- 
vironment modes. We are interested in the asymptotic 
behavior of this channel. That implies to take the limit 
n — > oo in the equations of Sec. III. It can be treated for 
some relations as the limit of Riemann sums resulting to 
the integral expressions. Instead of a set of equations 
on eigenvalues we get a set of equations on functions 
which are spectral densities for the involved (infinite- 
dimensional) matrices. Below we denote the spectral 
densities by the same symbols as proper eigenvalues, but 
written in calligraphic and replacing the mode number h 
by a continuous parameter £, i.e., i u h — * I u £, o q h — * O q £, 
etc. 

Suppose that all modes belong to the third stage, which 
holds true if (see Eq. (fig)) ) 

1 7 ? (27V + l) + (l- ?7 )(2^ onv + l)/o(25) . 1 

w := — — m ; — ; > 1, 

2|a| 7/ + (l-ry)(2iVe nv + l) 

(39) 

where Iq is the modified Bessel function of the first kind 
and zero-order. The lower bound C_ in this case is given 
by Eq. (|3"T|) after a formal replacement cosh s — > Iq (2s) in 
Eq. (f3"2")) . This example explicitly shows the possibility 
of an enhancment of the lower bound on the classical 
capacity with increasing degree of memory s. 

It is convenient to use the parameter £ as arising from 
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FIG. 5: Schematic representation of the "quantum water 
filling" for the model Qij = Sij+i + <5<,j-i- The angle £ 
parametrizing the spectral density corresponds to polar angle. 
White, grey and black sectors correspond to the first, second 
and third stages, respectively. Arrows show change of stages 
with increasing of N. The parameter r marks the points of 
stages transition. 



5 
3 
1 



:i/4 3t/2 37C/4 % 




1.3 
1.2 
1.1 

1 



/ 

/ 2 



ji/4 7t/2 3JC/4 \ 



71/4 ji/2 331/4 ^ 




3C/4 3C/2 37C/4 | 



the spectrum of V^ nv -matrix 



£«{ = ( Nenv + ^ I ' 



±2s cos £ 



(40) 



labeling both modes (if £ G [0,7r]) and eigenvalues (if 
£ e [0, 27r]). Plus and minus in Eq. (|40|) stand for it = g 
and u = p, respectively. Due to the mirror symmetry of 
eigenvalues (|40|) over quadratures, the symplectic spec- 
trum and the stages distribution have to be symmetric 
with respect to the point 7r/2, therefore we restrict our- 
selves to consider spectral densities only defined in the 
interval [0, 7t/2]. 

If w < 1 we can have one of the following stages distri- 
butions according to the results of numerical calculations 
for finite n (accordingly to the procedure of Sec. Ill): 

i) a mixture of second and third stages (2,3,2); 

ii) a mixture of second and first stages (2,1,2); 

Hi) all modes belonging to the second stage (2,2,2) 
which happens for a single value N 2 of the param- 
eter N, given s, r/ and N env . 

If N > N 2 or N < N 2 we have the (2,3,2) or (2,1,2) 
case with the center of the interval [0, ir] filled by the 
third or the first stage, respectively. We label the point 
of stages transition by r G [0,7r/2]. The possible stages 
distributions and dependence of r from iV are sketched 
in Fig|5J The continuity of the spectral density C U £ at 
points of stages transition r requires to hold A UT = O ut 
which can be rewritten as 



Here u = q gives X qr = 



x(t) = rfL UT + (\-T])Z UT . (41) 

_ c 2 S cost/ 2 ( see Eq (JX5J) for 
(2,3,2) and u = p gives l pT = 1/2 for (2,1,2) (we use 
different quadratures in these cases because of either C q £ 
or C p £ has a transition point belonging to the interval 
ofr). 



FIG. 6: Going from top-left clockwise the spectral densities 
V£,v^,C q £,T q £ (for Qjj = 5i,j+i + 8ij-i) are plotted vs the 
parameter £ for N = 0, 0.05, 0.67, 1, 2, 3.5, 6, 9, 11 (from bot- 
tom to top curve for quantities u^,C q ^,T q ^, and from top to 
bottom curve for quantity v^). Solid, dotted and dashed parts 
of curves correspond to third, second and first stages, respec- 
tively. Dash-dotted curve corresponds to the case of all modes 
belonging to the second stage. The values of other parameters 
used are N cnv = s = 1, 77 = 0.5. 



Then, the transcendent equation for x (see Eqs. (p?o]) 
and l|27p) can be rewritten as an equation for r 



V 



7T 7T 





+ i-,/ 


I" JO 


/o 



Epid^ = —x, 

7T 



(42) 

where (ti,t 2 ) is equal to (t, r) for (2,1,2) and to 
(71-/2, 7T - r) for (2,3,2). Moreover, x is given by Eq. ([41]) 
and T q £ is the spectral density for the second stage which 
can be found as solution of functional equation obtained 
from Eq. ([22]) or its approximations (see Eqs. ([28]) . ([29]) ) 
after the replacements discussed at the begining of this 
section. By substituting r = 7r/2 in Eq. (|4l?|) we find N 2 . 
Comparing it with the actual energy restriction we get 
correct stages distribution. Then, solving Eq. (|4"!?|) with 
the found stages distribution we arrive at r and x. Fi- 
nally, C_ is expressed through these parameters as follows 
(see Eqs. © and ©): 



C = 1 - ~T 3 
7T 



- ff ((l - T^env) 



9 ( V^«?- 5 



-9[ VO^Oi 6 -- 



where 



O qt = 772^ + (1 - 
°i* = fZrf + C 1 " 
r 3 is equal to 7r/2 for (2,1,2) and to r for (2,3,2). 



(43) 
(44) 
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FIG. 7: Exact solution, first-order and zeroth-order approx- 
imations for spectral density T q % vs £ for N = 1 (left) 
and iV = 0.01 (right). The values of other parameters are 
A'cnv = 0.5, s — 2.5, n = 0.95. Solid, dotted and dashed 
parts of curves correspond to third, second and first stages, 
respectively. Functions with maximum and minimum varia- 
tions correspond to exact solution and zeroth-order approxi- 
mation, respectively. 




FIG. 8: On the left, the quantity C_ is plotted vs s for values of 
r\ going from 0.1 (bottom curve) to 0.9 (top curve) with step 
0.1. The values of the other parameters are N — N anv = 1. 
Solid parts of curves correspond to the third and first stages, 
respectively. Dotted part of curves correspond to the second 
stage. On the right, the maximum of C_ over V em is plotted 
vs Menv for values of n — 0.1,0.5,0.9 going from bottom to 
top curve. Solid and dotted curves corresponds to Q — I and 



+ respectively. The value of the other 



parameter is N — 1. 



The "quantum water filling" effect for the considered 
model is shown in Fig[6] for symplectic spectral densi- 
ties T>£, and spectral densities C q £, T q ^. Graphs of 
Zg£ calculated through exact solution, zeroth-order and 
first-order approximations are shown in Fig[7J Despite 
some visible difference between exact and approximate 
spectral densities the corresponding symplectic spectral 
densities are almost equal, thus resulting to the differ- 
ence less than 0.05% between C calculated exactly and 
approximately. The small value of this difference comes 
from the fact that the Holevo-x has zero derivative over 
eigenvalues of V\ n and V c \ in the neighborhood of solutions 
of KKT-conditions (as they are equations for optimiza- 



tion problem). 

Finally, let us discuss the role of squeezing and mem- 
ory in lossy bosonic channel. Suppose, that the average 
amount of photons in the environment Af Gnv is fixed and 
we want to compare the quantity C for Q ^ I with that 
of O = I. As far as the Holevo-x is symmetric over 
modes one can expect that the lower bound on one shot 
capacity for O — I will always be higher. However, this 
is not true as results from the symmetry breaking over 
modes. Actually, we can see this from Figj8]-right where 
C_ maximized over parameters s, N env for memory and for 
memoryless cases is plotted versus M onv . In FigJ5]-left the 
lower bound on the classical capacity for memory model 
is plotted versus s for different values of rj. The limit 
of the one shot capacity (|57|) when s — > oo is still valid. 
Presumably, as we can see from Fig[8]-right, this limit 
is valid in the more general case of function maxy cnv (C) 
when M rriv — > oo. 



VI. CONCLUSION 

Summarizing, we have developed a powerful and ver- 
satile tool for the estimation of Gaussian quantum chan- 
nels capacities. It is based on the perturbative expan- 
sion of the von Neumann entropy of Gaussian states as 
function of the symplectic eigenvalues of the quadratures 
covariance matrix. We have applied this method to the 
classical capacity of a lossy bosonic channel. The found 
lower bounds are reliable rates of communication which 
can be achieved using solely Gaussian states for encoding 
classical information. By providing upper bounds it was 
shown that, in some settings, Gaussian state encoding is 
indeed optimal. Thus we have generalized the results of 
Q and partially those of Q and [l6j |. The expansion of 
the von Neumann entropy can be as well applied for the 
evaluation and optimization of other entropic function. 
Of a special interest is the case of the quantum capacity. 
In particular, it has been proven [?J that Gaussian state 
encoding is optimal for reaching the quantum capacity of 
the so-called degradable Gaussian channels. 
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